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Abstract — A geoacoustic inversion technique for high-frequency 
(12 kHz) multibeam sonar data is presented as a means to clas- 
sify the seafloor sediment in shallow water (40-300 m). The inver- 
sion makes use of backscattered data at a variety of grazing angles 
to estimate mean grain size. The need for sediment type and the 
large amounts of multibeam data being collected with the Naval 
Oceanographic Office's Simrad EM 121A systems, have fostered 
the development of algorithms to process the EM 121A acoustic 
backscatter into maps of sediment type. The APL-UW (Applied 
Physics Laboratory at the University of Washington) backscatter- 
ing model is used with simulated annealing to invert for six geoa- 
coustic parameters. For the inversion, three of the parameters are 
constrained according to empirical correlations with mean grain 
size, which is introduced as an unconstrained parameter. The 
four unconstrained (free) parameters are mean grain size, sedi- 
ment volume interaction, and two seafloor roughness parameters. 
Acoustic sediment classification is performed in the Onslow Bay 
region off the coast of North Carolina using data from the 12kHz 
Simrad EM 121A multibeam sonar system. Raw hydrophone data 
is beamformed into 122 beams with a 120-degree swath on the 
ocean floor, and backscattering strengths are calculated for each 
beam and for each ping. Ground truth consists of 68 grab samples 
in the immediate vicinity of the sonar survey, which have been an- 
alyzed for mean grain size. Mean grain size from the inversion 
shows 90% agreement with the ground truth and may be a useful 
tool for high-frequency acoustic sediment classification in shallow 
water. 





TABLE I 






Sediment Types 




Phi Value <j) 


Mean Grain Diameter 


Sediment Type 




(mm) 




<(-1.0) 


> 2.0 


gravel/rock 


(-1.0) -4.0 


0.06-2.0 


sand 


4.0-8.0 


0.004-0.06 


silt 


> 8.0 


< 0.004 


clay 



can range from clay (w 0.0039 mm) to boulders (w 256 mm) or 
greater. A phi value <fi scale conveniently represents the mean 
grain size according to 



log 2 



(1) 



I. Introduction 

The U. S. Navy has great interest in seafloor characteriza- 
tion due to its importance in shallow-water operations, such as 
landing operations, mine burial, and safety of navigation. De- 
termining a suitable route for communications cables, requires 
detailed knowledge of the seafloor and is another application 
for characterization of the ocean bottom. 

Obtaining and analyzing physical core samples or grab sam- 
ples provides an accurate characterization of the seafloor, how- 
ever, it is a time-consuming process and is not generally per- 
formed with sufficient coverage on an ocean survey. As an alter- 
native, acoustic seafloor characterization allows adequate cov- 
erage in much less time and, since sonar data is often collected 
on surveys, no additional data collection is required. The acous- 
tic data evaluated in this paper was collected in Onslow Bay 
with the 12 kHz Simrad EM 121 A Multibeam Echo Sounder. 

A. Sediment Types 

One of the most useful descriptors for bottom characteriza- 
tion is sediment type based on the mean grain diameter, which 



where d is the mean grain diameter in mm and do is the refer- 
ence diameter 1 mm. Approximate <fi values for selected sed- 
iments are given in Table U according to the Wentworth scale 

ID. 

B. Onslow Bay 

Onslow Bay off the coast of North Carolina is a challenging 
region for high-frequency acoustic sediment classification be- 
cause the bottom is dynamic (sediment drift) 1 2 1, heterogeneous 
in areas [3 1, |4| with shells, etc., mixed with the sediment, and 
is often composed of a hard bottom |4], J5J covered with only 
a thin (few centimeters or less) layer of sediment. The sonar 
data set from this survey is raw hydrophone data along the three 
parallel shiptracks depicted in Fig.^that are more or less par- 
allel to the coastline. Shiptrack 1 has 250 - 300 m water depth 
and is farthest from shore near the continental-shelf break. The 
seafloor slopes up to 0.5° in a direction perpendicular to the 
ship's heading. Shiptracks 2 and 3 are in shallower water (40 
- 60 m) about 80km from shore, seafloor here on the shelf is 
relatively flat. 

C. Hydrophone Data Processing 

Each port and starboard arrays are comprised of 64 hy- 
drophones. Each array is steered between -60 and 60 degrees 
(negative being in the port direction) in one-degree increments. 
However, both port and starboard arrays are steered at degrees 
so there are 122 steer directions (beams). 
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- Shiptrack 1, slope 

- Shiptrack 2, shelf 
■ Shiptrack 3, shelf 

□ Grab Samples 



-77.0 -76.8 -76.6 

Longitude (degrees) 

Fig. 1 

Shiptracks and collocated grab samples in the Onslow Bay 

AREA. 



1) Preprocessing: For each ping, header, raw data, and 
PAM (Power Amplifier Monitor) records are read from the tape 
or a file. A number of samples (usually 1028) from each ping 
are taken from the raw data record. Where applicable the be- 
ginning sample is selected according to the value of the Time- 
Varying Gain TVG and a hard-coded threshold. The effects 
of Programmable Gain PG, Fixed Gain FG, and Time-Varying 
Gain TVG are removed from the data. These computations are 
made in linear space based on values obtained from the header 
record. The data, now in units of digital number DN, are con- 
verted to sound pressure level SPL. Values for this conversion 
are taken from the header record. 

At this point the data are still at baseband. To beamform, 
the data are shifted to the original center frequency (12 kHz). 
To avoid aliasing the basebanded data must be resampled at 
a higher rate than the original sampling rate of 2.5 kHz. The 
resample factor used is 16, so the resampling rate is 40 kHz. 
The interpolation is done via a Fast Fourier Transform (FFT). 
The slow data are transformed to the frequency domain with a 
large FFT, shifted, then transformed back to the time domain. 

After the shift to 12 kHz, the average roll, pitch, heave and 
yaw for the given ping are computed. These values are then 
used to adjust the absolute locations (in software) of the receiver 
staves in the array and will enable the beams to be steered to 
consistent beam angles relative to the seafloor. 

Following the motion correction the data are beamformed by 
phase adjusting the frequency domain data according to the re- 
ceiver locations and desired steering angles. Taking the inverse 



FFT of these data yields a sound pressure P time series for each 
steering angle. The travel time of the bottom return is identified 
for each angle and an acoustic ray is traced out (here a constant 
sound speed profile is used because of a negligible sound speed 
gradient) to the corresponding bottom returns in order to obtain 
grazing angle. 

The sound pressure for the jth time sample of the ith beam 
is denoted P. L j . The data are converted to dB re /iPa and, based 
on the known geometry, the sonar equation is solved for bottom 
backscatter. 

2) Backscattering Strength: Backscattering strength BS is 
defined as 



BS = 10 log 



10 



(2) 



where is the backscattered sound intensity from an area of 
1 m 2 and / JIlc is the incident intensity at 1 m from the source 
|6|. The backscattering strength can be determined from the 
data by using the sonar equation 



BS = RL- SL + 2TL - I A, 



(3) 



where RL is the reverberation level (from the beamformed time 
series), SL is the source level, TL is the transmission loss in 
dB, and I A is the insonified area in dB re m 2 . The insoni- 
fied area is the area contributing to the received intensity and is 
computed using the 3 dB beam footprint, 



I A = 10 log^ J 2R 2 sin 



cot 



—cot 



or using the pulse length, 



.' ctR sin „ 
I A = 10 log] ' 



310 



cos 9 



smt> 
(4) 

(5) 



whichever is smaller, where R is the slant range to the bottom, 
8t is the transmit beam width, 9 r is the receive beam width, 
c is the water sound speed in m/s, and t is the pulse duration 
in s. The insonified area for several pressure time samples Py 
normally fall within the beam footprint, and the reverberation 
level for the ith beam RLi is averaged over these time samples, 



RU = 10 log! 



E 



Jo l J 



(6) 



Ji - h + 1 

where jo an d ji are the first and last time samples whose in- 
sonified areas lie within the ith beam's footprint. 

II. Backscatter Model 

The APL-UW backscatter model presented by Mourad and 
Jackson |7|, |8| treats the seafloor as a statistically homoge- 
neous fluid and predicts backscattering strength BS as a func- 
tion of grazing angle 6. The roughness of the bottom is de- 
scribed in this model by the bottom height spectrum. 
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TABLE II 



Model Input Parameters 


Parameter 


Symbol 


Description 


Density Ratio 


P 


density in sediment 
density in water 


Sound Speed Ratio 


V 


sound speed in sediment 

c 


Loss Parameter 


6 


imaginary wavenumber in sediment 
real wavenumber 


Spectral Strength 




Bottom height spectrum strength 


Spectral Exponent 


7 


Bottom height spectrum exponent 


Volume Parameter 




sediment attenuation coefficient 



(7) 



TABLE III 

Parameters Ranges 



where h is the reference height 1cm. The Mourad- Jackson 
model is valid for all frequencies between 10 and 100 kHz 
and is used here to represent the acoustic backscatter from the 
seafioor. 

Table [H] lists the six model input parameters, which, along 
with the sonar frequency / and sound speed c in water at 
the seafioor, determine both the roughness backscattering cross 
section a r (9) and volume backscattering cross section a v {9). 
The six input parameters are dimensionless except for (3 which 
has units of cm 4 . Combining these backscatter contributions 
from roughness (acoustic reflections from a randomly rough 
surface) and volume interaction (scattering of penetrating sound 
from sediment inhomogeneities) results in 



BS{9) = 101og 10 ((7 r + (7„). 



III. Data Inversion 



(8) 



The inversion problem is finding the set of input parameters 
that best fits the given data set. That is, which set of parameters 
minimizes the difference between the BS vs. 9 curve and the 
measured backscatter data. The sum of the squares of the data 
deviations from the model prediction is used as the measure for 
goodness of fit. 

A. Parameter Constraints 

If the six input parameters are unconstrained, the parameter 
space to be searched is six-dimensional. However, since cor- 
relations exist among some of the parameters, many solution 
parameter sets represent solutions that are physically unlikely. 
Hamilton and Bachmann |9|, |10| describe a relationship be- 
tween the parameters p and v and relate both to the mean grain 
size {(fi) of the seafioor sediments. Mourad and Jackson [7 1 pa- 
rameterize p, v, and 5 according to <f> values emphasizing the 
top few tens of centimeters of sediment, and the parameteriza- 
tion has been generalized to include coarse sand |8|. (Some 



Parameter 



Range 



<j> (-1.0) -9.0 

P within factor of 2 of APL-UW 

parameterization 1 8 1 
7 2.4-3.9 
a 0.00-0.02 



correlation exists between 5 and <f>, and the effect of physically 
meaningful values of S on the BS vs. 9 curve is negligible.) 
Gott 1 1 1 1 has used the idea of constraining some of the model 
parameters with some success. In addition the parameters used 
should be restricted to values that are physically likely. The 
parameter ranges used here are presented in Table IIIII 

The parameter space to be searched is now 4-dimensional (<j), 
(3, 7, a), and, since the backscatter model is highly nonlinear, 
one must be careful not to simply find one of the many local 
solutions. Two of the most common global search methods are 
simulated annealing and genetic algorithms. Both are suitable 
for most nonlinear problems. Simulated annealing (SA) is the 
best-fit search routine used here (e.g. see 1121 ). 

B. Simulated Annealing 

With the SA approach one searches the parameter space by 
continuously stepping to a new point in parameter space and 
computing the sum of the squares E for the data point residu- 
als. E is also known as the cost function. If the cost decreases 
from the previous location, the step is accepted. If, however, 
the cost increases, the step is only occasionally accepted. The 
probability p that a higher-cost step is accepted depends both 
on the amount of increased cost AE and on a variable referred 
to as temperature t according to the Boltzmann distribution, 

p = e- A£ /'. (9) 

This process is known as the Metropolis algorithm 1 13]. Local 
minima are escaped because of the steps of increased cost. The 
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temperature variable is gradually decreased until the probability 
of a higher-cost step is zero. The stepsize is also reduced slowly 
as the algorithm settles into the global minimum. 

IV. Results 

A. Slope Region 

The data for shiptrack 1 (farthest from shore) was grouped 
into bins of 200 pings covering an area of seafioor approxi- 
mately 3 kmxl km. The backscattering strengths in each bin 
were averaged according to grazing angle, and a best-fit pa- 
rameter set was found via simulated annealing for the averaged 
data. To illustrate Fig.[2]shows backscatter data for the first 200 
pings for shiptrack 1 along with the SA best-fit model curve. 
A comparison of inversion phi values with the analyzed grab 
samples is shown in Fig. [3] All 58 inversions for the 200-ping 
bins result in phi values indicating medium or fine grades of 
sand. The inversion phi values are in most cases only slightly 
greater than medium sand measured at the nearest grab sample 
location. 

B. Shelf Region 

Because of a higher ping rate, the backscatter data from ship- 
tracks 2 and 3 are binned in groups of 500 pings each. Figures^ 
and|3(closest to shore) compare the phi values from the SA in- 
version to the grab samples in the shallower water (m 40-60 m) 
with 6 1 of the 7 1 inversion phi values matching the nearest grab 
sample sediment type. Forty-three grab samples in the shelf re- 
gion show a medium or coarse sand bottom with one sample 
indicating gravel. The inversion yields 60 sand, 9 gravel, and 2 
clay values. This region exhibits greater variation in phi values 
than the slope region. 
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Fig. 4 

Comparison of phi value from inversion with grab sample 
analysis in the shelf region along shiptrack 2. 



V. Conclusions 

The inversion results are in good agreement with the ground 
truth (both in sediment type, i.e, sand, and in grade of sand) 
for the slope region where the sediment layer is known to be 
relatively deep and homogeneous (near the continental shelf 
break). The seafioor in the shelf area, on the other hand, is 
known to have little or no sediment layer and shells, rocks, etc., 
at the bottom. The inversion from the shelf region also agrees 
in sediment type with most grab samples, however, there is of- 
ten a discrepancy in grade of sand. Moreover, in a few cases 
the phi value from the inversion is the lower limit (-1, i.e., 
gravel/coarse sand). Because the sonar frequency is 12 kHz, 
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TABLE IV 

Average Phi Values 



Region 


Grab Samples 


Inversion 


Sediment Type 


% Agreement 


Slope 


1.92 ±0.36 


2.02 ±0.24 


Medium-fine sand 


100% 


Shelf 


0.82 ±0.71 


0.92 ±2.04 


Medium-coarse sand 


86% 



Clay 



7 - 

6 - Silt 

5 - 

4 

3 - 

2 - Sand 
1 - 

- 
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Fig. 5 

Comparison of phi value from inversion with grab sample 
analysis in the shelf region along shiptrack 3. 



the sound will penetrate any sediment layer less than about 13 
cm (wavelength) deep and interact with the hard subbottom. 
The backscattering strength predicted by the APL-UW model 
in this case will be invalid. 

Of the 129 inversions for the three shiptracks, 119 of them 
(92%) agree with the nearest grab sample in sediment type. Av- 
erage values and standard deviations are listed in Table llVl along 
with the percent agreement of the inversion phi values with the 
nearest grab sample. 

We believe the inversion method described here is promis- 
ing for determining sediment type in areas of relatively homo- 
geneous sediment and at least a few tens of centimeters deep. 
This process currently also provides an approximation for thin 
sediment layers or sediment with heterogeneous mixtures. 
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